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INSTANTONS IN THE QUANTUM FRAMEWORK OF 2D GRAVITY 
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We analyze the non— perturbative features of 2D quantum gravity defined by stochastic regularization of the unstable 
matrix model showing, first, that the WKB approximation of the well-defined quantum Fokker-Planck hamiltonian 
corresponds to the semiclassical eigenvalue density of the former. The double scaled potential exhibits an instanton- 
like behaviour, which is universal and scales, but whose interpretation in terms of pure gravity is still open. 



1) INTRODUCTION 

In the last year considerable effort has been 
devoted to solve the non- perturbative ambigu- 
ities in the possible definition of 2D quantum 
gravity trough matrix model techniques. While 
the solution[l] for the k — 3(fc odd in general) 
multicritical model is better understood [2], the 
problem for the k even models is much more com- 
plicated and it is still open whether the non- 
perturbative features of one-matrix models are 
related to pure gravity at all [3]. The difficulty 
is related to the instability of the matrix model 
with potentials unbounded from below, an old 
problem of euclidean gravity for which the two 
classical cures have been tried also here. The con- 
tour deformation [4] leads to a complex solution, 
presumably related to the complex Borcl sum 
of the perturbative expansion, whose imaginary 
part is proportional to the non-perturbative am- 
biguity obtained from the linearized Painleve 
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equation and has been interpreted as an instan- 
ton effect [6]. The second method is stochastic 
regularization, where one converts the ill defined 
D = problem into a well defined D = 1 — like 
one, and so it is considered by many to be the 
only realistic hope for a consistent real solution, 
albeit suspect as a possible source of additional 
arbitrariness. So far, only numerical computa- 
tions have been attempted [7] , some of them re- 
lated to the supersymmetry breaking [12] , which 
is a very interesting aspect of the original pro- 
posal [5] directly related to the D = 1 properties 
which will not be discussed here. 

Most of the numerical results rely on the de- 
coupled N-fermion formulation, which turns out 
to be a property of potentials of degree < 4. 
We prove here also that a mean field may be in- 
troduced in the WKB approximation to decou- 
ple the system, so that the semiclassical limit 
of the Matrix Model is reproduced. This way, 
the D — 1 system becomes an ideal Fermi gas 
whose spectral density is expressed in terms of 
the Fokker-Planck potential. 
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We arc then able to investigate analytically 
the leading non-perturbative behaviour of the 
stochastic stabilized quantum gravity also in the 
WKB approximation. We find that the double 
scaled (universal) stabilized Fokkcr Planck po- 
tential exhibits in the k — even case instanton- 
like effects, related to the non-Borel summa- 
bility and corresponding to the lifetime of a 
metastable state which is universal and char- 
acteristic of the non-perturbative contributions. 
On the other hand there are no such effects in 
the k = 3 multicritical potentials, as one should 
expect. 

On the light of this result we comment the pre- 
vious numerical attempts and discuss the leading 
effect of this tunnelling [bounce) contribution of 
the D = 1 — like spectrum in the free energy 
of the D = model, which can have different 
interpretations. 

2) WKB APPROXIMATION AND THE SEMI- 
CLASSICAL MATRIX MODEL 

Let us consider the D = hermitean one 
matrix model defined by the partition function 



D<f>e 



-pv(4>) 



(1) 



where <fi is an N x N hermitean matrix, and V(</>) 
is a generic potential of degree L: 



V{ct>) = Y,9nTr{ct> n ) 



(2) 



n=2 



This model can be understood in terms of the 
ground state of a D = 1 quantum mechanical 
system, whose hamiltonian is the (positive semi- 
definite) Fokker-Planck hamiltonian [11,8] 



H FP = Tr(P 2 ) + W FP 



d 

' d(j>ji 



W FP = 



1 dV dV (3 d 2 V 



4 d(j)ij d<j)ji 2 dfajdfyj 



(3) 



This means that the expectation value of an op- 
erator Q(4>) in the matrix model can be defined 
as the vacuum-expectation-value (VEV) of the 
quantum operator Q{<j)) 



(Q) = (0\Q(4>)\0) = J D(j)^l(<j))Q((j)) 
H FP ^ (4>) = E $ ((f)) 



(4) 



where E and ^o(0) are the energy and the 
wave function of the ground state. E$ = corre- 
sponds to the case when the potential is bounded 
from below, the matrix model is well defined, 
and *o(^) = exp(-(3V((f))/2); while eq.4, with 
E > 0, defines (Q) in terms of the true ground 
state wave function. 

Following well known techniques [9] , it is nat- 
ural to make a change of variables to the eigen- 
values Xi of the matrix <j>, by introducing the ef- 
fective ground state wave function 



*S // ({Ai}) = II( A *-^)*»W 

i<3 



(5) 



is totally antisymmetric, and it describes a 
gas of N Fermi particles. In general, the Fokker- 
Planck potential, Wfp, can be splitted into 
its diagonal (D) and non-diagonal (ND) parts, 
which after diagonalization of the matrix read 



FP 



w, 



(D) 



AND) 



w, 



(D) 



FP 



'FP +W FP 
pp. N 

xE((nA,)) 2 



(6) 



i 



-4X + n 9n^ 



n-1 



N 



n— 3 

L n-4 



n — 3 — s 
3 



i,j=l \n=4 



s=0 



Where X = N/(3 = e 70 is related to the (bare) 
2D cosmological constant in the usual way Ob- 
viously the ND piece does not vanish in general 
if L > 4, and the Fermi gas is not decoupled. 
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Nevertheless, in the semiclassical WKB approx- 
imation (/3 w fi," 1 — > oo), a mean field approx- 
imation (a la Hartree-Fock) may be performed 
to decouple the system. We show below that 

7>(0 fc )Tr(<^) « AT (uJkTr{(j) p ) + u> p Tr((f> k )— 

-Nio k io p ) + ■■■ (7) 

where the normalization is fixed by 

(Tr((b k )Tr(r))c « A^ 2 + 0(1/JV)) ; 

w fc = i < Tr(^ fc ) > (8) 

is consistent with the semiclassical limit of the 
matrix model. 

Under eq.7, the quantum mechanical system 
becomes an ideal Fermi gas of N particles, with 
the hamiltonian 

N 

H FP «]T/i FP (A0 

hFP W = -^ + tMA) 

U F p(X) 



8X 2 + 4 

I L-2 I L 

AX E E 39M-i-2 I A 



(9) 



i=l \.?'=«+2 



1 L \ 2 

+32 ~ o E E + (V'(X) f 

i=4 j=2 / 

Notice that the potential t/^p is bounded from 
below, and the one-fermion hamiltonian has 
a well-defined discrete spectrum, h FP (j) n (X) = 
e n <f> n {X). All the relevant information about the 
D = 1 quantum mechanical model is contained 
in the density, p(X, e) = (A|<5 (/i FP — e) |A), and 
its normalization fixes the Fermi energy, ep, 

= I dX I de p(X,e) (10) 



N 



The energy integral of p(X, e) provides the quan- 
tum mechanical version of the semiclassical den- 



sity of eigenvalues in the matrix model 

1 f ep N f 

MA) = -y de p(X, e) : X — — = J dX u(X) 

which, in the WKB approximation, reads 

U WKB {X) = _L^ eF _ UFp{X) e{ep _ Ufp{X)) 

Ztt 

1 (Tr Q(4>)) =JdX Q(X) u WKB {X) 



(11) 



Let us go back to the D = matrix model. We 
shall use the Schwinger-Dyson loop equations 
[10] to get a detailed expression for the semiclas- 
sical density of eigenvalues, and compare with 
the WKB approximation of the quantum me- 
chanical system, eq.ll. In the semiclassical (pla- 
nar) limit, the generating function of monomial 
expectation values, 

F{p) = -{Tr-^—) c (12) 
satisfies the Schwinger-Dyson equation 
F(pf - V'(p)F(p) + 

L-2 / L \ 

+ X E E isrt-i-j = » (13) 

»=0 \j=i+2 / 

where ujk are the "constants of integration" of 
this loop equation, and have already been defined 
in eq.8. Therefore, eq.13 is solved as 

,SC (A) 
X 



F(p) ee / dX^-t 
J P- 

= \ (V'( P ) y/Atfj) 



(14) 



L-2 / L 

A(p) = (V'(p)f IX Y, E 39i<*i-i-2 I P 

i=0 \j=i+2 

The "constants of integration", Wfe, are fixed by 
the condition that the imaginary part of F(p) 
defines a proper semiclassical density of eigen- 
values, u sc (X). Obviously, A(p) is a polynomial 
in p and all the branch cuts of eq.14 will be 
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squared root branch cuts. Therefore, A(p) has to 
satisfy the following constraints [9,4] : (i) A(p) 
must have only real zeros in the complex p-plane, 
and (ii) A(p) cannot have three consequtive odd 
degree zeroes. Under these conditions, the semi- 
classical density of eigenvalues is given by 

u sc (X) = -Im F{\) 

= 9(-A(\)) (15) 

The square root branch cuts of F(A), i.e., the in- 
tervals between odd degree (real) zeros of A(A), 
are the "bands" on which u sc (X) has support. 

Therefore, under the above mentioned restric- 
tions, the comparison between eqs.(14), (15) and 
(9), (11) shows that the WKB limit of the D = 1 
Fokkcr-Planck hamiltonian, with the mean field 
approximation of eq.7, reproduces the semiclassi- 
cal limit of the matrix model. Besides, the Fermi 
energy can be identified in terms of the uj k 

e L 

^ = .92 + 3 53 wi + (16) 

L I 1 i-2 \ 

+ Y1 + 2 51 (uj-iUi-j-i) i9i 

i=4 \ j=2 ) 

This result particularizes for the D = her- 
mitean one-matrix model the general arguments 
about the stabilization of bottomless euclidean 
field theories [13]. It also agrees with, and gen- 
eralizes, previous results obtained for the sim- 
plest potentials with L = 3 [5,11,8] and L = 4 
[12], showing that the critical behaviour of the 
Fokkcr-Planck hamiltonian is precisely that of 
the matrix model. 

Therefore, it is possible to describe the her- 
mitean matrix model like an ideal Fermi gas of 
N particles, whose one-fermion potential is fixed 
by the semiclassical density of eigenvalues 

U FP (X) =e F - [2iru sc (\)] 2 (17) 



This relationship formally holds only when the 
above mentioned constraints on A = U fp — &f 
are satisfied, and u sc is well defined. Never- 
theless, the quantum mechanical system, and 
u WKB , is defined even when this is not the 
case. Such quantum mechanical configurations 
arise when the naive ground state wave func- 
tion, = exp(— PV/2), is not normalizable in 
the WKB approximation. Therefore, they are re- 
lated to the true ground state with E > 0. 

3) CONTRIBUTION OF METASTABLE 
STATES AND INSTANTONS 

Once we have shown that the semiclassical 
density of eigenvalues is just the WKB density of 
the Fokker-Planck hamiltonian, we can use the 
matrix model results to get information about 
the Fokker-Planck potential Ufp- In particular, 
we want to find the features of Ufp that are 
directly related to the critical behaviour and, 
therefore, which will survive after the double 
scaling limit. 

The main problem of the matrix model for 
pure gravity is that the susceptibility which 
should be determined by the Painleve equation 
is only given as an asymptotic series, which is 
non-Borcl summable and so the difference be- 
tween any Borel sum (corresponding to different 
boundary conditions) and the hypothetical exact 
solution is proportional to T s e s . The 
exponential factor is an instanton-like effect as 
expected from the large order behaviour of the 
series, and it has been analyzed in the semiclas- 
sical limit of the matrix model and interpreted 
as the imaginary part of the complex Borel sum 
[6], related to the obstruction of real solutions 
[4], with the conclusion that this matrix model 
does not define the sum over topologies[15]. 

The first attempts to understand this prob- 
lem in the stabilizing Fokker-Planck [7] was the 
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computation of the D = free energy, which was 
obtained differentiating with respect to a source 
term added to the Fokker-Planck potential and 
nothing related to the instability was observed 
in the numerical results. This motivates further 
the present attempt to investigate these non- 
perturbative aspects analytically. 

Let us start by computing (F(p)) c , eqs.12, 14, 
using orthogonal polynomials [14] (we restrict 
ourselves to the case of even potentials): 

(F(p))c = (18) 



JV-1 



N-l oc 



' n=0 P H> yM n=0 j=0 1 



N-l 



= -y 



i 



dx- 



^ ^ - 4R(n) Jo yV - 4R(x) 

Taking the imaginary part, we get the semiclas- 
sical density of eigenvalues [9] 



uscip) 



Im{F(p)) c 



1 



x 9(4R(x)-p 2 ) 
dx- 



7T Wo jm(x) - p 2 

Next, we make the change of variables, x = 
W(R(x)), with W(R) defined in the usual way 
in terms of the potential 

an the final expression reads 



usc(p) 



1 



R(X) W'(R) 

dR [ ' = 6{AR{X) -p 2 ) 



71" Jp*/i y/AR - p 2 



^~V e F - U F p(p)0(e F - U F p(p)) 

Z7T 



(19) 



As a first result, we see that the semiclassical 
range for the eigenvalues is given by p 2 < 4R(X). 



Now we can obtain the critical behaviour of 
uscip) m the double scaling limit. For the fcmo- 
del it corresponds to #1 



W{R) = X c (1 - 7 (i? c - R) k ) 
X = X c {l- /r^T) 
R = i2c(l /(T)) 



(20) 



Therefore, 



V4fl(X)-p2 



usc(p) = ^- k X c Jdy(AR c -p 2 -v 2 ) k - 1 (21) 

Notice that, from eqs.19, 21, U F p(X) has a root 
of order 2k — 1 at A 2 = 4R C , in the critical point, 
X = X c and R = R c . 

Let us specialize eq.21 for k = 2, pure gravity. 
Then, 



uscip) = ~rX c (AR c -p 2 )^m(X)-p 2 - 

47T L 



--(AR(X)-p 2 )i 



(22) 



This means that the critical behaviour of the 
Fokker-Planck potential is 



U FP (X) - e F = 



(23) 



1 X C 



1 X C 



(A 2 - 4R(X)) (6R C - 2R(X) - \ 2 f 



(A 2 - m c f 



120- 



7 



-(A 2 — AR C )T 



Now, the leading behaviour of the potential 
around Ao = ±\/AR c when — > oo is 

^ P (A)«±16r«7^(^-^) (24) 

# x Notice that X = and, with these general defin- 

tions, 

1 fT\ i 



f(T) 



R< 
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where z = [3$ (A — A ). Therefore, the double 
scaling limit of Ufp has two (symmetric) degen- 
erated secondary minima at z^ = ± y/T/ (AjRc), 
if T > 0, and one absolute minimum at A = 0. 

The perturbative expansion around A = 0, the 
absolute minimum of the potential, corresponds 
to the semiclassical WKB expansion, which, as 
it is well known, reproduces the \/N expan- 
sion of the matrix model[13]. In this expan- 
sion, the subdominant terms in eq.24 are cru- 
cial because the normalization condition of hav- 
ing N — > oo bound states depends also on (3 
(in fact, these terms ensure that the potential 
has bound states!). Nevertheless, there is also 
a non-vanishing contribution related to the sec- 
ondary minima which implies the existence of a 
metastable state, corresponding precisely to the 
Fermi level, which decays due to barrier pene- 
tration, and whose non-perturbative lifetime re- 
mains finite in the double scaling limit. Such life- 
time is propotional to the inverse of the imagi- 
nary part of the metastable state energy, as we 
show below. 

The dominant barrier penetration contribu- 
tion can be computed in the semiclassical ap- 
proximation [15]. We consider the hamiltonian 

h FP {\) = --^ + Plu FP (\). The correspond- 
ing euclidean partition function can be written 
as the L — > oo limit of the functional integral 

tr (e- LH ) = f [dq(t)]e- s ^ (25) 

Jq(-L/2)=q(+L/2) 

where S[q(t)] is the euclidean action 

r+L/2 / 1 p,2 \ 

S[q(t)}= J ^ dt[-?(t) + ^U(q(t))) 

and tr{e~ LH ) is given by exp—LE n , where 
the sum extends over the whole spectrum. There- 
fore, in the L — ► oo limit, 

Im tr(e~ LH ) oc Im e~ L{E ™ + iE '^ 

cx Le- LEr -El (26) 



where he have taken into account the smallness 
of the, non-perturbative, imaginary part of the 
energy E l m of the lowest eigenvalue that becomes 
complex. 

We assume that the potential has a secondary 
minimum at q = q m , which obviously requires 
that the equation U(q) = U{z m ) has, at least, 
another solution; we shall call qo the one that 
is closest to the minimum, assuming qo < q m . 
Therefore, the would-be eigenvalue correspond- 
ing to the secondary minimum becomes com- 
plex, providing the dominant contribution to 
Im tr(e~ LH ). Its imaginary part is given by an 
instanton, which is a solution of the euclidean 
equations of motion that starts from the sec- 
ondary minimum at euclidean time L = — oo, is 
reflected in qo, and comes back to the minimum 
at time L = +oo. 

The euclidean equation of motion is 

\qc=^U'{q c ) (27) 
which has a conserved quantity 
q 2 c - P 2 U{q c ) = -B 2 U{q rn ) = -E r m (28) 
Therefore, the instanton action is 

r +L/2 /i b 2 \ 

%J = J dt^-q 2 (t) + ^-U(q c )j 

r+L/2 1 

= LE r m + dt-q 2 c - LE r m + (29) 

J — L /2 

+ dqy/p{U{q) - U{q m )) = LE r m + S t 
J qo 

The integration around the saddle point corre- 
sponding to q c in the gaussian approximation 
provides the instanton contribution to the par- 
tition function. This calculation involves a Jaco- 
bian factor, which becomes complex because of 
the turning point q , and is proportional to L 
because of traslation invariancc in the euclidean 
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time. The final result is 

Im tr(e- LH ) oc Le~ s ^ = Le~ LE " l ~^ i 

oc Im exp(-L(E r m + ie~ Si )) (30) 

Therefore, e~ Si corresponds to the imaginary 
part of the eigenvalue whose real part is the value 
of the potential at q m , E = E r m + ie~ Si , which 
corresponds to a metastable state 

Let us particularize this formula to the case of 
pure gravity, eqs.(9,24), where 

U(q) =f3-iv(z) ; v(z) = Az 3 - BTz 

q = ^/W c + {}-iz (31) 

and E r m is just the semiclassical value of the 
Fermi energy. Then, Si becomes finite in the dou- 
ble scaling limit (fi — > oo) and the result is 

Therefore, the lifetime of this metastable state is 

Toce-** ; S^^X^RD^Ti (33) 

As mentioned above, no related effect has been 
observed in the numerical computation [7] . This 
negative result could be related to the numerical 
precission . In [7] the the D = specific heat is 
computed by introducing a source in the Fokker- 
Planck potential, Ufp(^) + J\ 2 . In this way, the 
derivative of the corresponding free energy with 
respect to the source provides < Tr $ 2 >, which 
is related to the first derivative of the pure grav- 
ity free energy, < Tr $ 2 >oc F' = - < P >. 
The consistency of this method requires that the 
value of J is small enough to keep the phase 
structure of pure gravity. Implementing this for 
the generic even matrix model potential of fourth 
degree shows [16] that in the double scaling limit, 
[3 — ► oo, this requires that the source scales, 

4 

J = (3~sj and if this scaling is not explicitly 
performed, the constraint on the source J fixes 



the precission to be a function of T decreasing 
substantially at small values of T. In particu- 
lar, for the range explored in [7], to preserve the 
precision achieved at large T when T is small, 
one should increase the value of N with a factor 
of three. This is precisely the region where the 
instanton-like behaviour should be observed. 

Let us finally consider eq.21 for the next mul- 
ticritical model, k = 3. Notice that in general, 
eq.21 only gives the leading terms of use in the 
double scaling limit, but it is exact when the 
canonical representative is chosen for the poten- 
tial. When k = 3 this means a potential of degree 
6. In such case, the result is 

u sc (p) oc (4R(X) - p 2 y ((p 2 + R(X) - 5R C ) 2 + 
+5 (R(X) - i? c ) 2 ) (34) 

which means that, in this case, Ufp(^) has only 
one absolute minimum at A = 0. Therefore, 
no metastable states are present in the k = 3 
model, and no imaginary part is expected. This, 
of course, agrees with the fact that the k = 3 
model is already well defined as a matrix model. 

We have been considering non -pcrturbative 
effects in the D = 1 — like spectrum of the dou- 
bled scaled Fokker-Planck potential which sta- 
bilizes the matrix model, whose original Borel- 
nonsummability is refleceted in the appearance 
of real instantons from the secondary minima 
of the potential. The problem is how to trans- 
late these results into the actual D = prob- 
lem. This requires the computation of the non- 
perturbative contributions to the D = 1 density, 
not yet completed, which leaves the interpreta- 
tion of this instability open at present. On the 
other hand, one can clearly stablish [16] that the 
tunneling scales properly, remaining finite in the 
limit and that it involves the instanton action in 
eq. 33. 
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4)C0NCLUSI0NS 

We have analyzed the problem of the non 
perturbative definition of 2D pure gravity in k— 
even matrix models with the Fokkcr-Planck sta- 
bilization, discussing other proposals. We have 
proved that the fermionic formulation, on which 
previous results were based, is only valid for 
k < 4 , but that with a mean field in the WKB 
approximation the generic potential decouples 
and it is explicitly related to the density of states 
of the matrix model. We have shown then how to 
analyze analyticallly the critical behaviour of the 
corresponding Fokker-Planck potential perform- 
ing explicitly the double scaling limit It turns 
out that the scaled potential has secondary min- 
ima in the relevant k even case, which exhibit 
instanton-like behaviour corresponding to the 
Fermi level becoming complex and which could 
reflect the non-Borel sumability of the pure grav- 
ity series solution of the Painleve equation. We 
have discussed how this result can be seen in one 
of the previous numerical computations and why 
is not seen in the other. The final interpretation 
of these instabilities in terms of sum over sur- 
faces, or the pure gravity true vacuum, requires 
more work as mentioned at the end of last sec- 
tion, which is in progress. 
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